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Abstract 

A metastable state, or a false vacuum, is not an eigenstate of the Hamiltonian in quantum 
field theory. Its energy density has a non-zero imaginary part equal to its decay width. 
Therefore, supersymmetry cannot be exact in the false vacuum. We calculate the size of this 
effect using the path integral approach. 
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The Universe is believed to have arrived at its present state via a series of phase transi- 
tions, each of which involved a decay of some metastable quantum state. A good understand- 
ing of these transitions is essential for explaining the physical world at present. In particular, 
the baryon asymmetry of the Universe, dark matter, and many other cosmological and ex- 
perimental data may owe their existence to the physics of the early Universe. Although the 
experimental evidence for supersymmetry has yet to be discovered, theories with low-energy 
supersymmetry provide an attractive explanation for the hierarchy of scales and are widely 
considered viable. It is of interest, therefore, to understand the effects of supersymmetry on 
the decay of a metastable state. 

It is well-known that supersymmetry is broken at non-zero temperature M. There is, 
however, an additional effect, a generic feature of the false vacuum which is due to its 
instability per se, which can contribute to SUSY breaking^] and be equally or more important 
than the finite-temperature effects. 

Unlike the true vacuum, |0), a false vacuum |0') is not an eigenstate of the Hamiltonian 
H [0]. It's decay width per unit volume, (T/V), is equal to the imaginary part of the 
energy density, Im(0'\H\0'), which must be non-zero if the false vacuum is to decay. On 
the other hand, by the usual argument, exact supersymmetry cannot be realized unless the 
vacuum energy is zero. The latter requirement implies that both the real and imaginary 
parts of the energy density must vanish. Indeed, the anticommutation relation between the 
supersymmetry generators Q a 

{Q a ,Q\} = 2a^P lx , (1) 
together with the statement of unbroken supersymmetry, 

Q a \0') = Ql\0') = 0, (2) 

application to the false vacuum, the term SUSY breaking is somewhat a misnomer because super- 
symmetry is not really broken. The unstable state that is considered to be a "vacuum" is simply not a part 
of the spectrum and the corresponding Fock space is ill-defined. Nevertheless, for all intents and purposes 
SUSY will appear broken to an observer in the false vacuum, which justifies the use of the term. 
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yield immediately that 

(0'|P M |0') = /m(0'|P M |0') = o r/v = 0, (3) 

in contradiction with the assumption that |0') is a metastable state and thus has a non-zero 
decay width. 

This simple line of reasoning shows that supersymmetry cannot be exact in any metastable 
state and that the vacuum in which SUSY is unbroken must be stable. The latter can serve 
as an independent argument for the stability of the supergravity vacuum with vanishing 
cosmological constant even in the presence of lower lying minima in the potential, a fact 
established in Ref. || by other means. 

It is clear that the size of the SUSY breaking effects in the false vacuum is of order 
r/V. Below we will show that this is, in fact, the case and will compute the mass splitting 
Am 2 = m 2 — m 2 ~ (F) for the bosonic and fermionic components of a chiral superfield 
which arises due to non-perturbative effects related to false vacuum instability. In general, 
there are other contributions to Am 2 which can be computed in perturbation theory. For 
example, if the scalar potential is of the form Fig. |I|, and if supersymmetry is unbroken in 
the true minimum, then there will be a contribution (Am 2 ) R to Am 2 from the real part of 
the energy density due to spontaneous or explicit (soft) SUSY breaking at <fi = 0. In realistic 
models one can assume (Am 2 ) R ~ M 2 [/sy . At finite temperature supersymmetry is always 
broken []T| which gives rise to a temperature-dependent contribution (Am 2 ) T . So, in general, 

(F) ~ Am 2 = (Am 2 ) s + (Am 2 ) T + (Am 2 ), (4) 

where (Am 2 ), is the non-perturbative piece induced by the false vacuum instability. The 
relative magnitudes of the terms in (|j) are model-dependent. However, one can always 
construct a model in which the barrier separating the false vacuum from the true is sufficiently 
small, so that the last term in equation (^) is arbitrarily large, while the first two terms are 
much smaller. 



2 



We see that the supersymmetry breaking due to the instability of the false vacuum will 
dominate over that due to finite-temperature and other effects for a metastable state that 
is sufficiently short-lived, or, in other words, if T/V ^> and T/V ^> T 4 . In what 

follows we will neglect the (Am 2 ) R and (Am 2 ) T which can be calculated perturbatively in 
each given model and will concentrate on (Am 2 ) r 

We now evaluate the mass splitting of the fermion tp and boson components in a model 
involving a single chiral superfield $ = {<j),ijj,F} with the scalar potential that has two 
non-degenerate local minima (Fig. [1]). 

In the path integral formulation a vacuum expectation value of some function f(<p) of 
the dynamical field can be written as 

(/(«/>)> = lim^-r In / [d4] exp(-~S[0] + J Jf (</>)), (5) 

where S[4>] is the Euclidean action and J(x) is the source. 

The non-perturbative contribution to /(</») can be evaluated in the semiclassical limit 
h — > using the saddle point approximation (see, e.g., [0, U). The idea of the method is 
that one can expand the path integral around the finite-action solutions which are saddle 
points of the Euclidean action and are, therefore, determined by the variational equation 

5S = 0. (6) 

For a fermion ip(x), the Euler-Lagrange equation (Euclidean equation of motion) cor- 
responding to equation ([|) is a first-order differential equation. The boundary conditions 
ip(x — > oo) = are necessary to ensure the finiteness of the action. The only solution ip{x) 
of equation (|]) satisfying these boundary conditions is the trivial one ip(x) = 0. 

In contrast to fermions, the bosonic equation of motion corresponding to (|^) is second- 
order and it allows for a non-trivial solution of finite action, the so called "bounce" [0. 
The bounce dominates the path integral contribution to the imaginary part of the energy 
density, which was calculated in Refs. 0, [|, More specifically, the bounce for a scalar 
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u(4>) 




Figure 1: The false and the true vacua. 

potential U(<p) (Fig. 1) is a nontrivial 0(4)-symmetric |J field configuration determined by 
the equation 2 



A0(r) = -tu(^) 



with the following boundary conditions: 



{ (d/dr)0(r)\ r=o = O 
6(oo) = 



where r = a/x 2 . 

The potential U(<f>) can be written as 



Z n>2 



and the boson mass, by definition, is 



(7) 



(9) 



2 If the potential depends on several scalar fields, then finding the bounce numerically becomes a non- 
trivial task. In this case, the method of Ref. Q is useful. 
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je- s ^{d(j)\ 

= m " + }eW/ <10) 
where f((f>) = (d 2 /d<p 2 ) £ n>2 c„0 n . We note that /(0) = 0. 

We now apply the saddle point method to calculate the quantity Am 2 = — rriQ. First 
we will compute the normalization factor 

N = f e" SM [#] (11) 



in the denominator of the second term in (|T0|) . This part of the calculation will follow 



closely that of Refs. [Q, [|]. For simplicity, we consider the case of one space-time dimension 
t; generalization to d = 4 is straightforward and will follow. 

The bounce of equation (|7]) is centered at t — because (fi'(0) = 0. Let's consider a 
contribution from a single bounce (fi(t — tj) centered at t = tj, — T/2 < tj < T/2, where T is 
some large quantity taken eventually to infinity. A classical field configuration corresponding 
to n widely separated bounces will also be a saddle point of the integrand with the Euclidean 
action nS[(j)]. Therefore, we have to sum the contributions of all such saddle points and 
integrate over the positions of the bounces. 

We are interested in the limit T — ► oo. The n bounces centered at t\, t^-, ■■■ t n , where 
—T/2 < t n < ... < t 2 < t\ < T/2, are separated by vast regions where the field is nearly 
zero. Therefore the determinant part of the path integral can be written as a product 



CO 



/uh)e- wT ' 2 K n (12) 



where u = U"(0) and K is defined so as to yield the correct answer for n = 1. 
We must now integrate over the locations of the centers of the bounces. 
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T/2 rti rtn-i T n 

dti / dt 2 ... / dt n = — (13) 

-T/2 J -T/2 J -T/2 TV. 



The summation of contributions with different n gives: 



OO _ rp n 

(4)1/2 exp( _^ T/2 + Ke- s ^' h T) (14) 

This is the same as the corresponding result in 0]. The quantity K is the correction to the 
energy which determines the transition probability. We will come back to it later. 

Now we treat the numerator of the second term in fllPp. This is basically the same, except 
now the configuration which comprises n widely separated bounces, = S?=i <K^o ~tj)i wm 
contribute 

n 

^/[^-t^e-W (15) 

where the property /(0) = is essential. 

Like before, we have to integrate over the positions of the bounces: 

rT/2 rh rtn-i JL _ 

/ dh / dt 2 ... / dt^f^to-tj)] (16) 

J -T/2 J-T/2 J -T/2 ~{ 

To evaluate the integral, we exchange the order of summation and integration. Then we 
let each variable ti run over the entire interval [—T/2, T/2]. The latter change would lead 
to a double counting of the contributions differing by permutations of the identical bounces. 
We therefore divide by (n!), the number of permutations of n bounces: 



I dt i I dt ? ■■■ I dt n 4>(t -tj) 

r i n 



,T/2 


fT/2 


fT/2 


/ dh 


/ dt 2 ... 




J-T/2 


'-T/2 


J-T/2 
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1 n r T/2 r T/2 _ r T/2 

*i .- { / *i/[«<b - «,)]} / *. (17) 



J'=l' 



In the last expression, the integration over i,- gives (in the limit of large T) /f^ </>(£)d£, while 
each of the other integrals gives a factor of T. The sum over j has n identical terms and is 
equal to n times the value of each. The factor n then cancels with that of n\. Altogether we 
get 

/•oo _ f>n-1 



We notice in passing that the expression (18) differs from (13) in two respects: (1) there is 



a constant factor of / f[(j)(£)]d£; (2) it has one less power of T. 

And, finally, we sum over n: The n = component will not contribute; therefore, the 
sum starts from the n = 1 term: 



i j C OO - '/'it — J. 

E(^) 1/2 { / fmme- T ' 2 K n {e- s ^Y^— 

n=l ™ J-oo [n ~ 1) 

/OO — / i 
-oo 7m , 



oo _ rpn—L 

u \l/2 „-uT/2 T sn-l( „-S[4>]/%\n-l 1 



if f Ke- S[ ^ h {^-fl 2 exp(-ujT/2 + Ke- sm/h T) (19) 

J-oo Tin 

As we now divide expression fllTf ) by fll4j) , we obtain a simple result: 



/OO _ 
/[0(O]<£} ^e" 5M/?i (20) 
-oo 

where the factor K is that computed in |4]][]: 

3 In applications, one is often interested in transitions between phases with different symmetries. In 
this case, some additional zero modes appear in the determinant of equation (|2l]) due to the Goldstone 
phenomenon. The corresponding N zero eigenvalues should then be omitted from the "primed" determinant 
in (21) and an additional factor [J ff{x)d i x\ N ^ 2 will appear in the expression for Am 2 ||. 



7 



iAm 2 i = {j nm m ^ e-^ 



Sm\-sum det'[-dl + U»$)] " J 2 



x(i + o(n)) (2i) 



det[-#j + [/'< 

where det' stands for the determinant from which all zero eigenvalues are omitted. 

As was explained earlier, there is no analog of the bounce for fermions, so that the 
correction fl2~T| ) applies to bosons only. 

The above discussion dealt with the zero-temperature tunneling. It is straightforward to 



generalize the expression ([H]) to the finite-temperature case ||. The changes will amount to 
replacing the Euclidean action of the bounce with the action of a three-dimensional bounce 
63 divided by the temperature, and the factor (S/h) 2 with (S^/h) 3 / 2 /T 1 / 2 . 

In summary, we have computed the mass splitting between the components of a chiral 
superfield which is related to the instability of the false vacuum. If the metastable state 
decays sufficiently fast, this effect may dominate over the SUSY breaking finite-temperature 
effects, as well as the soft SUSY breaking terms of phenomenologically acceptable size. 
We have also given an independent argument for the supersymmetric vacuum stability in 
super gravity. 
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